Rounding of First Order Transitions in Low-Dimensional Quantum Systems with 

Quenched Disorder 
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We prove that the addition of an arbitrarily small random perturbation to a quantum spin system 
rounds a first order phase transition in the conjugate order parameter in d < 2 dimensions, or for 
cases involving the breaking of a continuous symmetry in d < 4. This establishes rigorously for 
quantum systems the existence of the Imry-Ma phenomenon which for classical systems was proven 
by Aizenman and Wehr. 
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A first order phase transition, in Ehrcnfcst's terminol- 
ogy, is one associated with a discontinuity in the density 
of an extensive quantity. In thermodynamic terms this 
corresponds to a discontinuity in the derivative of the 
free energy with respect to one of the parameters in the 
Hamiltonian, more specifically the one conjugate to the 
order parameter, e.g. the magnetic field in a ferromag- 
netic spin system. In what is known as the Imry-Ma phe- 
nomenon 

[Ell, 

any such discontinuity is rounded off in 
low dimensions when the Hamiltonian of a homogeneous 
system is modified through the incorporation of an ar- 
bitrarily weak random term, corresponding to quenched 
local disorder, in the field conjugate to the order param- 
eter. 

This phenomenon has been rigorously established for 
classical systems [H, E| , where it occurs in dimensions d < 
2, and d < 4 when the discontinuity is associated with 
the breaking of a continuous symmetry. In this letter 
we prove analogous results for quantum systems at both 
positive and zero temperatures (ground states). 

The existence of this effect was first argued for random 
fields by Imry and Ma on the basis of a heuristic analysis 
of free energy fluctuations. While the sufficiency of Imry 
and Ma's reasoning was called into question, the pre- 
dicted phenomenon was established rigorously through a 
number of works 0, 0, H, IS 0, @, 0| • The statement was 
further extended to different disorder types by Hui and 

id, El- 
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The general existence of the Imry-Ma phenomenon in 
quantum systems was not addressed by these rigorous 
analysis, and in particular the Aizenman- Wehr 0, 0| 
proof of the rounding effect applies only for classical sys- 
tems. However, as stressed in [ljj], establishing whether 
the Imry-Ma phenomenon extends to first order quantum 
phase transitions (QPTi) is an important open problem. 
The results presented here answer this question. We find 



that the critical dimensions for the phenomenon for quan- 
tum systems are the same as for classical systems, includ- 
ing at zero temperature. 

We consider spin systems on the d-dimensional lattice 
Z d , where the configuration at each site is described by 
a finite-dimensional Hilbert space, with a Hamiltonian of 
the form 



(1) 



where {k x } are translates of some local operator Kq, and 
h and e are real parameters. The quenched disorder is 
represented by {t] x }, a family of independent, identically 
distributed random variables. TLq may be translation in- 
variant and nonrandom, or it can include additional ran- 
dom terms (although we will not discuss the latter case, 
our results hold there also). For convenience we will as- 
sume that ||/t x || = Ij which can be arranged by rescaling 
h and e. We will refer to the rjs as random fields, al- 
though in general they may also be associated with some 
other parameters, e.g. random bond strengths. 

An example of a system of this type (with k x = a x ) 
is the ferromagnetic transverse-field Ising model with a 
random longitudinal field [3] (henceforth QRFIM), with 
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where <j x ^ (i — 1, 2, 3) are single-site Pauli matrices, and 
Jx-y > 0- The QRFIM has recently been studied as a 
model for the behavior of LiHo y Y i X F4 with x > 0.5 in 
a strong transverse magnetic field [lj, LL5J . 

We will examine phase transitions where the order pa- 
rameter is the volume average of the expectation value of 
k x with respect to an equilibrium state (that is, a state 
satisfying the Kubo-Martin-Schwinger condition 
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and show that this quantity cannot be discontinuous in 
h for low-dimensional systems. As is well known, this 
order parameter is related to the directional derivatives 
(±) of the free energy density, 

m±(T,h,e) := -g^C^M) (3) 
where, as usual, at positive temperatures 

?(T,h,e)= lim —LlogTre-^ r (4) 

ryz d p\L | 

(with (3 := 1/ksT), and ^"(0, h, e) is the corresponding 
limit of the ground state energy. Here Hr is the Hamilto- 
nian of the system restricted to the finite box T C 7L d , and 
|T| is the number of sites in that box. It is known under 
the assumptions enumerated below that for almost all r\ 
this limit exists and is given b y a non-random function 
of the parameters (see, e.g. 0, 17|), which does not de- 
pend on the boundary conditions. By general arguments 
which are valid for both classical and quantum systems, 
T is convex in h; therefore the directional derivatives ex- 
ist, and are equal for all but countably many values of h 

For typical realizations of the random field, the in- 
terval [m-(T,h,e),m + (T,h,e)] provides the asymptotic 
range of values of the order parameter for any sequence 
of finite volume Gibbs states or ground states (the argu- 
ment is similar to that found in [4| for classical systems). 
At a first order phase transition m_ < m + , and there 
are then at least two distinct infinite volume equilibrium 
states [16[ with different values of the order parameter. 
In the QRFIM the m+ and the m_ states can be obtained 
through the + or — boundary conditions (i.e. the spins 
a x 3 ^ are replaced by ±1 for all x £ T). In general, such 
states are obtained by adding ±5 to the uniform field h 
and letting 6 — > after taking the infinite volume limit. 

Our discussion is restricted to systems satisfying: 

A. The interactions are short range, in the sense that 
for any finite box A £ 1 d the Hamiltonian may be 
decomposed as: H = Ha + Va + Ha? , with Ha acting 
only in A, Ha" only in the complement A c , and Va of 
norm bounded by the size of the boundary: 



||Va|| < C\dA\ . 



(5) 



B. The variables r\ x have an absolutely continuous distri- 
bution with respect to the Lebesgue measure (i.e. one 
with a probability density with no delta functions), 
and a finite rth moment, for some r > 2. 

Our main results are summarized in the following two 
statements. The first applies regardless of whether the 
order parameter is related to any symmetry breaking. 

Theorem 1. In dimensions d < 2, any system of 
the form of ([TJ satisfying the above assumptions has 
m+(T, h, e) = m_(T, h, e) for all h, and T > 0, provided 



The next result is formulated for situations where the 
first order phase transition would represent continuous 
symmetry breaking. An example is the O(N) model with 



Ho = - 22 Jx- 



y u x 



(6) 



where a are the usual quantum spin operators. More 
generally, 7io is assumed to be a sum of finite range terms 
which are invariant under the global action of the rotation 
group SO(N), and a x is a collection of operators of norm 
one which transform as the components of a vector under 
rotations. With the random terms the Hamiltonian is 



H = Ho-J2ft 



(7) 



Theorem 2. For the SO(N)- symmetric system de- 
scribed above, with the random fields ff x having a 
rotation-invariant distribution, the free energy is contin- 
uously differentiable in h at h = whenever e ^ 0, d < 4, 
and N > 2. 

Before describing the proof, let us comment on the 
implications of the statements, and their limitations. 

1. While the statements establish uniqueness of the 
expectation value of the bulk averages of the observables 
k x , or a x (in Theorem 2), they do not rule out the possi- 
bility of the coexistence of a number of equilibrium states, 
which differ from each other in some other way than the 
mean density of n, which they share. More can be said for 
models for which it is known by other means that non- 
uniquess of state is possible only if there is long range 
order in k. (Such is the case for QRFIM, through its re- 
lation to the classical ferromagnetic Ising model in d + 1 
dimensions [18||.) 

2. The results address only the discontinuity, or sym- 
metry breaking (as in the QRFIM), but they leave room 
for other phase transitions, or singular dependence on h. 
For instance, for the Ashkin- Teller spin chain for which 
Goswami et. al. [13] report finding the Imry-Ma phe- 
nomenon in some range of the parameters but not else- 
where, the results presented here rule out the persistence 
of a first-order transition between the paramagnetic and 
Baxter phases in the full range in the model's param- 
eters. However, they do not rule out the possibility of 
other phase transitions. 

3. Randomness which does not couple to the order 
parameter of the transition need not cause a rounding 
effect. For example, in the transverse- field Ising model 
in a random transverse field, where the random field r/ x 
couples to cr^, ferromagnetic ordering is known to per- 
sist , . Presumably the same is true for the Baxter 
phase of the Ashkin- Teller model. It was even suggested 
that there are systems in which the introduction of ran- 
domness of this sort may even induce long range order 
which would not otherwise be present 2^, 21 1, and our 
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results do not contradict this. In addition we can draw no 
conclusions about quasi-long-range order, that is power 
law decay of correlations, including of k x . 

Other comments, on the technical assumptions un- 
der which the statements hold, are found after the proofs. 

The proofs of Theorems [1] and [2] are based on the analy- 
sis of the differences, between the m + and the m_-states, 
in the free energy (at T = 0, ground state energy) which 
can be ascribed to the random field within a finite region 
A of diameter L. Putting momentarily aside the ques- 
tion of existence of limits, a relevant quantity could be 
provided by: 



lim lim Av \G S A vIvaI 



Av [G 



A,r] 



(8) 



where G S A r is the difference of free energies 



GiAv) 

with 
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,h+S 



-ir),h—5 



W, h-S 



ylogTr exp(~(3H^ h ), 



(9) 



(10) 



7/ A ) is the random field configuration obtained from r\ by 
setting it to zero within A, and Av [• |t?a] is a conditional 
expectation, i.e. an average over the fields outside of A. 
(The modification of the field h by ±<5 serves to select 
the desired (m±) states). 

Somewhat inconveniently, it is not obvious that for all 
models the limits in (|9|) exist. Nevertheless, one can prove 
that for each system of the class considered here there is 
a sequence of volumes Tj / Z d for which the limit exists 
for all A, with convergence uniform in r\ A . The proof 
of this assertion is by a compactness argument, whose 
details can be found elsewhere 



22] 



The essence of the proof of Theorem 1 is the contra- 
diction between two estimates: 



Under Assumption lAl equation ([S]): 
\G a (t])\ < 4C|9A| . 



(11) 



ii. Whenever m_ < m + , Ga/v|A| converges in distri- 
bution to a normal random variable with a positive 
variance (as one would guess by considering the dif- 
ference in the random field terms between states of 
different mean magnetizations, neglecting the states' 
local adjustments to the random fields). 

More explicitly, for the upper bound we note that in 
the absence of the interaction terms V A , the right hand 
side of ([9]) would be zero. Using (O, one gets (fTT|) . 

To prove the normal distribution for G Ar , we apply 
a theorem of 0, Proposition 6.1] (as corrected in 23, 



p. 124]). It implies that for A /* 1 d , under Assump- 
tion [B] G A /y/\K\ converges in distribution to a normal 
random variable with variance of the order of 



b = Av 



dG A 
df] x 



e(m+ — mJ) 



(12) 



The two statements described above contradict the 
assumption that ro_ < m+ in dimensions d < 2. That 
is so even at the critical dimension, where L d / 2 = L d ~ x . 
The reason is that the lower bound implies the existence 
of arbitrarily large fluctuations on that scale, whereas 
the upper bound is with a uniform constant. This proves 
Theorem 1. 

The above proof is similar to that of the classical 
results 0, 0] which this work extends. However, the 
discussion of the free energy fluctuations was based 
there on the analysis of the Gibbs states, and more 
specifically of the response to the fluctuating fields of the 
'metastates' which were specially constructed for that 
purpose. Except for special cases, such as the QRFIM, 
that argument was not available for quantum systems, 
where the equilibrium expectation values are no longer 
given by integrals over positive measures. The proof of 
the quantum case is enabled by a more direct analysis of 
the free energy. 

Theorem 2 is proven by establishing that in the pres- 
ence of continuous symmetry the upper bound (|11[) . for 
A = [— L, L] d , can be replaced by: 



\G A (r) A )\<KL 



d-2 



(13) 



Here G A is defined as in d9j) , dHJ , but h = 0, and S is 
replaced by by 8 := Se with e a unit vector. This change 
in the upper bound raises the critical dimension to d = 4. 
To obtain (11311 we focus on 



Av K 



? iJ,Se 



f), — 5e 



'/A 



(14) 



Since G A (rj A ) = lim^o hm r ^z<* \ (sa.i^a) ~ 3A,r( )) > 
any uniform bound on |<7 A r | for given A implies a similar 

bound on Ga|. The claimed bound may be obtained 
though a soft-mode deformation analysis, which we shall 
make explicit for the case of pair interaction (the general 
case can be treated by similar estimates). 

The free energy F^' e in Q14p may be rewritten by 
rotating both the spins and the field vectors with respect 
to an axis perpendicular to e at the slowly varying angles 



0, < L 

Mp^Tr, L<\\x\\< 2L 
7T, ||x|| > 2L 



(15) 



The rotation aligns the external fields in the two terms 
(±5e), except in the region ||x < 2L\\ where the effect 
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is negligible when S — > 0. The effect of the rotation on 
the random fields is absorbed by rotation invariance of 
the average. In the end, the Hamiltonian of the rotated 
system differs from the Hamiltonian used to define the 
other free energy by 



{x,y}CT 



(16) 



When the resulting expression for F ] 



ff,—Se 



in (U4J is 



expanded in powers of 9 X — 9 V 



IT IX 



y\\/L, the 



zeroth-order term cancels with F^'~ Se , and the second 
and higher order terms yield the claimed bound. The 
main difficulty is to eliminate the first order terms, which 
amount to a sum of 0(L d ) quantities each of order l/L. 
However, the sign of these terms is reversed when the 
rotation is in the reversed direction. To take advantage 
of this, we combine two expressions for g A p(^a) with the 
rotations applied in opposite directions, yielding: 



giAVA) = Av [logTre- 011 
-ilogTre-^+ A ^ 



ilogTr e 



-/3(ff+Aff_e) 



VA 
(17) 



(where H = H^ ,Se .) By known operator inequalities [l^ |: 
logTr e-^-ilogTr e -^ H+AH ^ - MogTr e -0(*+Afr_ 4 



< ±\\AHg + AH-e 



(18) 



The right hand side is zero to first order, and one is 
left with an upper bound on gj^ r of the desired form. 
Repeating this analysis with the roles of the terms 
exchanged we obtain an identical lower bound, and thus 
inequality (fT3|) follows. 



The above argument is spelled out in detail in [22|. 
Let us end with few additional comments on the assump- 
tions. 

4. For Thcorcm[2j the assumption that the interaction 
has a strictly finite range can be weakened to a condition 
somewhat similar to Assumption[A] For pair interactions 
(equation 0) it suffices to assume: 



\J X \ \\x\\ 2 < oo. 



5. The restriction to absolutely continuous distribu- 
tion excludes a number of models of interest. Such an 
assumption is generally necessary at zero temperature, 
as can be seen by the behavior of the Ising chain in a 
random field [24| which takes only a finite number of val- 
ues. For positive temperatures it can be replaced by the 



requirement that the distribution has a continuous part 
which extends along the entire range of values. For the 
QRFIM at finite temperature one need only assume that 
the random field has more than one possible value, and 
this may well be the case more generally. 
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